(N 

o 

(N 



C/2 



C^ 



- 1 - 

Streaming cold cosmic ray back-reaction and thermal 
instabilities 

Anatoly K. Nekrasov^ and Mohsen Shadmehri^ 

^ Institute of Physics of the Earth, Russian Academy of Sciences, 123995 Moscow, 
Russia; 



CN . anekrasov@ifz.ru, nekrasov. anatoly ©gmail. com 

"<C ! ^ Department of Physics, Golestan University, Basij Square, Gorgan, Iran; 

T^ ' m.shadmehri@gu.ac.ir 
O 



ABSTRACT 



We investigate the streaming and thermal instabilities of the electron-ion plasma with 
homogeneous cold cosmic rays drifting perpendicular to the background magnetic field 



> 

m 

[^^ ' in the multi-fluid approach. One- dimensional perturbations along the magnetic field are 

in ■ 

^ ' considered. The induced return current of the background plasma and back-reaction of 

o' 

CN ' cosmic rays are taken into account. It is shown that the cosmic ray back-reaction results 

in the streaming instability having considerably larger growth rates than that due to the 
rS • return current of the background plasma. This increase is by a factor of the square root 



of the ratio of the background plasma mass density to the cosmic ray mass density. The 
maximal growth rates and corresponding wave numbers are found. The thermal instability 
is shown to be not subject to the action of cosmic rays in the model under consideration. 
The dispersion relation for the thermal instability includes ion inertia. In the limit of fast 
thermal energy exchange between electrons and ions, the isobaric and isochoric growth rates 
are derived. The results obtained can be useful for the investigation of the electron-ion 
astrophysical objects such as galaxy clusters including the dynamics of streaming cosmic 
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rays. 
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1. INTRODUCTION 

There is a growing interest towards understanding of interactions of cosmic rays 
with plasma systems in astrophysics and their possible effects. Irrespective of various 
mechanisms which are proposed for the generation of such high energy particles, cosmic 
rays may interact with the existing turbulent motions in a plasma or even may excite 
them. In order to study cosmic rays, a particle description is needed, although the fluid 
approximation is used as well for simplicity. Cosmic rays may further induce ionization 
which may dramatically change the physical properties of a system. For example, the 
ionization by cosmic rays has a vital role in star formation near the Galactic center (e.g., 
Yusef-Zadeh et al. 2007) or in the dead zone of protoplanetary disks (Gammie 1996). On 
the other hand, the heating rate is enhanced because of cosmic rays and this important 
effect has been studied in the context of structure formation in the interstellar medium 
(ISM) via thermal instability (e.g.. Goldsmith et al. 1969; Field et al. 1969). 

Another contribution of cosmic rays to the dynamical evolution of the system is their 
pressure. Many authors studied the dynamical role of cosmic rays in structure formation at 
large scales by Parker instability (Parker 1966; Kuwabara & Ko 2006), magnetorotational 
instability (Khajenabi 2012) and even galactic winds or outflows (e.g., Everett et al. 2008). 
Recently, Wagner et al. (2005) and Shadmehri (2009) extended the classical thermal 
instability (Field 1965) to include cosmic rays. The thermal instability has been used to 
explain existence of structures not only in the ISM but also at the very large scales like 
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in the intracluster medium (ICM). This instabihty is assumed as a possible mechanism 
for producing molecular filaments (Sharma et al. 2010) seen in galaxy clusters with short 
(<1 Gyr) cooling times (e.g., Conselice et al. 2001; Salome et al. 2006; Cavagnolo et al. 
2008; O'Dea et al. 2008). Recent linear analysis of thermal instability with cosmic ray 
pressure shows that the instability is suppressed (e.g., Shadmehri 2009). Cosmic rays has 
been included by Shadmehri (2009) (see also Sharma et al. 2010) in the framework of 
the magnetohydrodynamic equations as a second fluid having the velocity of the thermal 
plasma. Numerical analysis shows that the cosmic ray pressure can play an important role 
in the dynamics of cold filaments making them much more elongated along the magnetic 
field lines than the Field length. This is consistent with observations as well (Sharma et 
al. 2010; see also Snodin et al. 2006). Also, including cosmic rays is required to explain 
the atomic and molecular lines observed in filaments in clusters of galaxies (Ferland et al. 
2009). 

However, there is another important effect of cosmic rays which has not been considered 
in the context of thermal instability, namely the presence of streaming cosmic rays. These 
particles are charged and their drift motion induces a current. It has been assumed that at 
some distance from the shock cosmic rays are mainly positively charged particles (Riquelme 
& Spitkovsky 2009). The cosmic ray drift driving a constant current results in arising of 
the return current provided by the background plasma (e.g.. Bell 2004, 2005; Riquelme & 
Spitkovsky 2009). The possible role of this effect in the generation of thermal instability 
needs to be considered. There is also another important issue like the amplification of 
magnetic fields. The classical cyclotron resonant instability has been proposed long time 
ago to explain this process (Kulsrud & Pearce 1969). However, this mechanism has turned 
out to be unable to provide sufficient energy in the shock upstream plasma. In order to 
resolve this problem, just recently a new non-resonant instability has been introduced that 
may provide a much higher energy (Bell 2004). This instability, which is known as the Bell 
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instability, has also been confirmed by non-linear numerical simulations (e.g., Riquelme & 
Spitkovsky 2009). Subsequent works extended this instability into various directions by 
considering various physical factors (e.g., Reville et al. 2007; Reville & Bell 2012). However, 
all these works, except the paper by Bell (2005), have been restricted to the cosmic ray drift 
speed parallel to the initial magnetic field and correspondingly to unmagnetized cosmic 
rays, where Larmor radii of cosmic rays defined by longitudinal velocities are much larger 
than the length scale of interest (see also Zweibel 2003). In the case of magnetized cosmic 
rays, whose Larmor radii are smaller than the typical length scales of the system, the cosmic 
ray current can be perpendicular to the initial magnetic field as has been demonstrated by 
Riquelme & Spitkovsky (2010). Riquelme & Spitkovsky (2010) studied this perpendicular 
current-driven instability in the linear regime and also numerically and compared their 
growth rate with the cosmic ray current-driven instability (CRCD) by Bell (2004). The 
growth rates and length scales in both cases were similar. But the authors have not included 
the cosmic ray back-reaction analytically. 

Depending on the magnetic field strength and the cosmic ray fiux, the cosmic ray 
thermal pressure effect may become important and modify the maximum growth rate of the 
CRCD (for a detailed discussion, see Zweibel & Everett 2010). The thermal instability in 
galaxy clusters in the multi-fluid approach has been considered by Nekrasov (2011, 2012). 
The related effects by cosmic rays were not included in these papers. Here, we take into 
account streaming cold cosmic rays. We consider a geometry in which homogeneous cosmic 
rays drift across the background magnetic fleld and perturbations arise along the latter. 
Such a geometry is analogous to that treated by Riquelme & Spitkovsky (2010). However, 
we include the cosmic ray back-reaction. We also take into account the plasma return 
current (e.g.. Bell 2004, 2005; Riquelme & Spitkovsky 2009, 2010). For simplicity, we ignore 
the action of gravity here (Sharma et al. 2010). The effects of the gravitational fleld have 
been investigated in detail in multi-fluid approach in papers by Nekrasov & Shadmehri 
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(2010, 2011). Thus, our present study extends previous analytical studies by considering 
not only the thermal effects but the currents driven by cosmic rays and their back-reaction. 

The paper is organized in the following manner. Section 2 contains the fundamental 
equations for plasma, cosmic rays, and electromagnetic fields used in this paper. Equilibrium 
state is discussed in Section 3. In Sections 4 and 5, the perturbed velocities of the ions 
and electrons and perturbed plasma current are given, respectively. Corresponding results 
obtained for cosmic rays are contained in Sections 6 and 7. The total perturbed current is 
given in Section 8. Wave equations are found in Section 9. Dispersion relation including 
the plasma return current, cosmic ray back-reaction, and terms describing the thermal 
instability is derived in Section 10. In Section 11, a discussion of important results obtained 
is provided. Possible astrophysical implications are given in Section 12. Conclusive remarks 
are summarized in Section 13. 



2. BASIC EQUATIONS FOR PLASMA AND COSMIC RAYS 



The fundamental equations for a plasma that we consider here are the following: 

^"^- + V, . Vv, = -^^ + ^E+^v, X B, (1) 



dt ''^j^j ^i ^c 



the equation of motion, 

_1 + V . n,v, = 0, (2) 

the continuity equation, 

(9T- 1 

--^ + V, • VT, + (7 - 1)T,V . V, = - (7 - 1) -A (n„T,) + ul {n,,T,) {T, - T,) (3) 
ot rii 

and 

--^ + Ve ■ VTe + (7 - 1) TeV ■ Ve = - (7 - 1) - V • q, - (7 - 1) -C, (n^, T,) (4) 

ot Tie He 
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are the temperature equations for ions and electrons. In Equations (1) and (2), the subscript 
j = i,e denotes the ions and electrons, respectively. Notations in Equations (l)-(4) are 
the following: Qj and rrij are the charge and mass of species j, Vj is the hydrodynamic 
velocity, rij is the number density, pj = UjTj is the thermal pressure, Tj is the temperature, 
z/fg(ne, Te) (z/gj {rii, Te)) is the frequency of thermal energy exchange between ions (electrons) 
and electrons (ions) being vl^{ne-,Tf.) = 2z/je, where z/jg is the collision frequency of ions with 
electrons (Braginskii 1965), n^z/fg (ne,Tg) = Uei^ei (^i?^e), 7 is the ratio of the specific heats, 
E and B are the electric and magnetic fields, and c is the speed of light in a vacuum. Here, 
for simplicity, we do not take into account collisions between the ions and electrons in the 
momentum equation. This effect for the thermal instability has been treated by Nekrasov 
(2011, 2012), where, in particular, conditions allowing to neglect the collisions have been 
found. However, the thermal exchange should be included because its time scale compares 
with the dynamical time. The value qg in Equation (4) is the electron heat flux (Braginskii 
1965). In a weakly collisional plasma, which is considered here, the electron Larmor radius 
is much smaller than the electron collisional mean free path. In this case, the electron 
thermal flux is mainly directed along the magnetic field, 

qg = -Xgb (b ■ V) Tg, (5) 

where Xe is the electron thermal conductivity coefficient and b = B/i? is the unit vector 
along the magnetic field. We only take into account the electron thermal flux given by 
Equation (5) because the longitudinal ion thermal conductivity is considerably smaller 
(Braginskii 1965). We also assume that the thermal flux in the equilibrium is absent. 
The cooling and heating of plasma species in Equations (3) and (4) are described by 
function Cj{nj,Tj) = njAj (Tj) — rijTj, where Aj and Tj are the cooling and heating 
functions, respectively. The form of this function is somewhat different from the usually 
used cooling-heating function £ (Field 1965). Both functions are connected to each other 
via the equality Cj {nj,Tj) = nijnjfj. Our choice is analogous to those of Begelman & 
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Zweibel (1994), Bogdanovic et al. (2009), Parrish et al. (2009). The function A^ (T,) can 
be found, for example, in Tozzi & Norman (2001). 

Equations for relativistic cosmic rays we take in the form (e.g. Lontano et al. 2002) 

' + Ve^ • V {RcrPcr) = + Qcr ^+-Vcr X B , (6) 



at rirr \ c 



^--)(¥' 



'^cr 



^ + v,..V)(^^^)=0, (7) 



where 



r T 

Iter — ^ I 



In these equations, Pcr = 7cr^crVcr is the momentum of a cosmic ray particle having the rest 
mass rricr and velocity v^r, Qcr is the charge, Pcr = Icr^ncrTcr is the kinetic pressure, ricr is the 
number density in the laboratory frame, Tcr is the adiabatic index, 7^^ = (1 — v^^/c^) 
is the relativistic factor. The continuity equation is the same as Equation (2) for j = cr. 
Equation (8) can be used for both cold nonrelativistic, T^r <^ iricrC^, and hot relativistic, 
Tcr ^ iTT-crC^, cosmic rays. In the first (second) case, we have F^ = 5/3 (4/3) (Lontano et 
al. 2002) . The general form of the value Rcr applying at any relations between Tcr and 
rricrC^, can be found e.g. in Toepfer (1971) and Dzhavakhishvili and Tsintsadze (1973). 

Equations (l)-(4), (6), and (7) are solved together with Maxwell's equations 

VXE.-If ,9, 



and 



An. IdE 



where j = j^^ + jcr = Xli Ij^j^j + jcr- Below, we consider at first an equilibrium state in 
which there is a stationary cosmic ray current. 



3. EQUILIBRIUM STATE 

We will consider a uniform plasma embedded in an uniform magnetic field Bq (the 
subscript here and below denotes background parameters) directed along the z-axis. We 
assume that the plasma is in equilibrium and penetrated by the uniform beam of cosmic 
rays having the uniform streaming velocity Ucr along the y-axis. It is reasonable to suppose 
that in such uniform model the magnetic field due to cosmic rays will be absent. This 
picture is analogous to the consideration of the gravitational potential in the equilibrium 
state in an infinite uniform medium having a constant mass density. Then we obtain from 
Equation (10) 

J2 li^jo^io + J^^o+^^ = 0. (11) 

j 

From Equation (1), we easily find in the equilibrium state 

Eq xBq , c gEp 

^eO = C -2 , Vio = VeO H ^— , VeOz = ViQ^ = 0, (12) 

where we have assumed that d/dt ^ Ucj, oJcj = qjBQ/mjC is the cyclotron frequency. 
Analogously, we obtain from Equation (6) under condition Rcr'jcrd/dt <C oJccr {^ccr = 
qcrBo/mcrc) 

VcrO = VeO + Ucr- (13) 

In Equation (13), we have neglected the polarizational drift of cosmic rays (the second 

term on the right hand-side for the ion velocity in Equation (12)). This can be done due to 

condition riio ^ ^cro (see below), which is always satisfied. Using Equations (12) and (13), 

we will find the current jo 

qiHiocdEo , 

Jo = 5 ^— + gcr^crOUcr, (14) 

where we have taking into account the condition of quasi-neutrality qiniQ + qeneo + qcrficro = 0. 
The number density Ucr is the one in the laboratory frame. We note that this condition 
is satisfied in astrophysical plasmas due to cosmic ray charge neutralization from the 
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background environment (Alfven 1939). Substituting Equation (14) into Equation (11), we 
obtain 

I ^-+1 1 -Kj- + 4:7rqcrncrOUcr = 0, 

where cai = (Bq / AnmiriiQ) is the ion Alfven velocity. Usually, the inequality c ^ c^j is 
satisfied. Thus, the induced electric field Eq is determined by equation 

—-2 = _47rg^^ncrOUcr^- (15) 

ot c^ 

We note that this equation has been given in the paper by Riquelme & Spitkovsky (2010). 
Substituting Equation (15) into Equation (12) for ions, we find the return plasma current 

jret = QiniO (Vjo — Vgo) = —QcrncrO'^cr = " jcrO; (16) 

which is equal to the cosmic ray current and has an opposite direction. The induced plasma 
current drift velocity Upi = Vjq — Veo is equal to Upi = — (gcr^cro/O'i^io) ^cr- Using Equation 
(15), we see that the polarizational cosmic ray drift velocity {Rcr'^crc/ oOcrBo) (9Eo/(9t can be 
neglected in comparison with Uc,. under condition rriiniQ ^ Rcrlcrf^cr-ncro- If cosmic rays 
are not too relativistic, this condition is satisfied. The drift velocity Upi will be taken into 
account at the consideration of thermal instability which is provided below. 

We will here consider the case in which background temperatures of the electrons and 
ions are equal between each other, Tgo = Tjq = Tq. However, for convenience to follow the 
symmetric contribution of the ions and electrons, we keep in general calculations different 
temperatures. In this case, the thermal equations in equilibrium are given by 

C-i {riio, Tio) = Ce (neo, Teo) = 0. (17) 
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4. PERTURBED VELOCITIES OF IONS AND ELECTRONS 

Here, we will investigate one-dimensional perturbations depending on the 2;-coordinate. 
Equations for perturbed velocities of ions and electrons are given in the Appendix A. We 
consider Equation (A3) under condition w^- ^ d"^ /dt^ . Then, the transverse velocities are 
given by 



gj F _q^dEi^ qj d^Eiy 

^j^y -'-'Ix ~r 9 oj. "•" ■^ Oj.9 ' 

where Equation (A2) has been used. We also have taken into account that Biz = (see 
Equation (9)). The longitudinal velocities Vnz and Veiz are given by Equation (A13) 

Lviiz = Ha, Lveiz = Hei, (19) 

where the values L, Hn, and iJei are defined by Equations (A9), (A12), and (A14)-(A20). 



5. PERTURBED PLASMA CURRENT 

It is known that contribution of streaming fiows to the dispersion relation is due to 
the difference of their velocities (e.g., Nekrasov 2008, 2009a, 2009b, 2009c). Therefore, we 
can not take into account the contribution of the electric drift VgQ. We also could choose 
the appropriate frame of reference and make the electric drift initially zero (Riquelme & 
Spitkovsky 2010). 

The linear components of the plasma current j^n = ^ • qjnjoVji+ ^ . qjTijiVjQ have the 
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form 



4vrj,/i. -[l^— E,y + -J—- - -^^7^, (20) 



^ a;e, ; ^^ uji, dt coi dt^ 



4:7ijpny --{ 2^-— ]Eix + -T^S^ + 



^—^ oJh / "^ wE (9t c<;£ (9t2 






AnqiUioUpi- — 



1 fdY^dHii 



5 



47rjp/i2 = 4:TTqinio^ + Airqeneo^, 

where Upj = [Annjoq'j/mj^ is the plasma frequency and Upi = Vioy. When deriving 
Equation (20), we have used Equations (18) and (19) and the continuity equation. 



6. PERTURBED VELOCITY OF COSMIC RAYS 

We assume cold, nonrelativistic, Tcr <C TricrC^, cosmic rays here for which Tcr = 5/3 
(Lontano et al. 2002). In this case, the value Rcr is equal to the unity, Rcr = 1 (see Equation 
(8)). The interaction of cosmic rays with the background plasma can be considered 
without taking into account the cosmic ray back-reaction and including the latter. If 
cosmic rays are stationary in both equilibrium and perturbations, then we have to take the 
quasi-neutrality condition in the background state in the form gjUjo + 5'e^eo = to avoid 
the appearance of the current due to the electric drift. When the cosmic ray back-reaction 
is involved, the condition of quasi-neutrality becomes gjUjo + Q'e^eo + qcrficro = 0. The last 
condition have been used in Section 3. We here consider the back-reaction of cosmic rays 
also in perturbations. We assume that the condition u"^^^ ^ 7cro*^^/^''^ i^ satisfied. For 
ultrarelativistic cold cosmic rays {'-fcr ^1), the last condition can be violated. Then the 
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velocities Vcrix,y given by Equation (B3) are the following: 



_ Qcr rp QcrlcrO dEjy QcrlcrQ ^ ^Ix 

Vcrly -l^lx H 5 TT: I ^ ?77o ; 

rricrUJccr ^cr^ccr "^^ ^cr^ccr "^^ 

where we have used Equation (B2). 

The longitudinal velocity Vcriz is defined by Equations (B6) and (B7) 

Licr'^crlz -t^crl- \^^) 

Substituting Equation (21) for Vcriy into Equation (B7) for Hcri and using Equation (B2), 
we obtain 

H^,^ = _Jf^c' 7^ !^_^ (_E^^ + jcr^dEiy\ ^qc^d_ / ^ _ ^^\ ^ .23) 

rricrCOccr '^^'^ '^'^^ C^ Ozdt \ ^ Uccr dt J Ulcr Ot \ ^ C ^ ) 



7. PERTURBED COSMIC RAY CURRENT 

The linear perturbations of the components of the cosmic ray current jcri = 
gcr'^croVcri + Qcrncri^cr are equal to 

^TTJ^rlx = -—Ely + — ^ -5 --^, (24) 

Uccr ^ccr ^^ ^ccr '^'^ 

. ^pcr r-i , ^pcrlcrO oEiy ^pcrlcrO 9 Ei^ 



1 fdY'dH„i 



- ATTqcrUcrOU 
^T^jcrlz = ^-nqcrricrO 



cr 



Lcr \dt J dz ^ 

Her I 



T ' 

^cr 

where Upcr = {Aimcroqlr / "^cr) is the cosmic ray plasma frequency. When obtaining 
Equation (24), we have used Equations (21) and (22) and the continuity equation for cosmic 
rays. The value Hcri is defined by Equation (23). 
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8. PERTURBED TOTAL CURRENT 

We find now the components of tlie perturbed total current ji = jp/i + jcri- Adding 
Equations (20) and (24), we obtain 

< . _ I ^pi ^pcr'ycro\ dEix ^ f ^pi ^pcrlcro\ 9 Ely 






1 f d\~^ dHii , 1 /d\~^dHcri 



- AnqiUioUpl— I — 1 — AnqcrTlcrOUcr T \af I n ' 

4vrji^ = 4:7rqinio—— + Anqeneo—— + A-KqcrUcro^ — , 
wfiere we liave used tlie condition of quasi-neutrality g^njo + qe^eo + Qcr-ncro = 0. 



9. WAVE EQUATIONS 

To obtain the wave equations, we have to substitute the current (25) into Maxwell's 
equation (10). Omitting small terms under condition d/uddt <^ 1 and also assuming that 
7f^o djujccrdt ^ 1, we find in the one-dimensional case 

c'[l^~'^^-E. = e..Eix. (26) 

'dt ) "dz^ ~ Ey = ^yy^iy " ^yzEu, (27) 

= —EzyEly + EzzElz. (28) 
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Here, the following notations are introduced: 



^pi ^pcricrO 

'XX 2" ' 2 ' 



(29) 



CJ, 



-yy 



-yz 



'zy 



pt (^pcrlcrO / 2 ,2 



UjliUpi 



2 
cc 

DGi 
L 



DG2 2 2 1 

~r ^ pernor j 



d\-^ d^ 



dt 



dz^ 



, 2 2. 

~r ^pcr'^cr j 



dt d^' 



2D ( qeUeO^ A , 2 1 



pi 



L 



t-pl 



per 



qiUio 

£zz = {^peG?, + WpjCl) — + tOp- h 1. 



d_ 

dz'' 



L 



U 



When deriving Equations (26)-(28), we have used Equations (9), (23), and (A18) without 
the contribution of terms proportional to Vqx- 

Equation (26) describes magnetosonic waves including the contribution of cosmic 
rays at conditions under consideration. Equations (27) and (28) describe the streaming 
instability due to the cosmic ray flow and thermal instability subjected to an influence of 
cosmic rays. We note that when Ucr = 0, the thermal instability is defined from Equation 

(28), e,, = and Ei, ^ (Nekrasov 2011, 2012). 



10. DISPERSION RELATIONS 



Equations (27) and (28) are given in their general form which permits us to investigate 
different particular cases. Making use of the Fourier analysis for perturbations proportional 
to exp {ikz — iut), we obtain 



-1 



u^ 



^yy^zz 



^yz^zy 



(30) 
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Substituting into Equation (30) expressions for Eyy, Eyz, Szy, and Ezz, which are defined by 
Equation (29), we find 

^zz 5 1 2" 77 = ^^^'^1 "~ '^2tt3, (31) 

where Ca is the Alfven velocity including the contribution of cosmic rays 



In Equation (31), we have introduced notations 



2 2 -^^2 2 2 / \ 

"1 = ^piUpl—j— + UJpcrUcr-j—: (32) 



,2 



DGi 



"2 = UJ^^Upl—— + Up^^Ucr-j—, 
Li Li(^,f. 

2 ^ ( n , qeUeO^ \ , 2 1 

tt3 — ^pi~rUpl I (-72 H Ct4 I + UJ„^j.Ucr-p — • 

To calculate the right-hand side of Equation (31) with values 0:1,2,3 given by Equation 
(32), it is convenient to consider the expression L^ {szz(^i — <^2Q;3)- Carrying out the 
calculations and taking into account that qiiiiQ + qeUeo ~ and Ucr ^ ViQy, we obtain 

L^ {szzOii - 0203) = ^liuli {ul^G2G3 + WpjG'iG'4) D^ (33) 

+ ^pcr'^cr i^pe^S + ^piGl) D — . 

It can be shown that the value ^"^^20^ + cjpjG'iG'4 acquires the simple form 

a;J,G2G3 + ujIG^G^ = coll, (34) 

where we have taken g, = —qe (see Equation (A19)). The value oJp^Gs + ojpjGi can be given 
in the following form: 

uj%G, + ujIG, = -u^D [u' - k'Gl) , (35) 

where C^ is given by 

(36) 
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In the value ezz on the left-hand side of Equation (31), we neglect the contribution of 
the unity which has arisen due to the displacement current. It is easy to see that it can 
be done if iJi^ ^ w^ for the cold plasma and bJi^ ^ ^^"^Te f*-*^ "^^ warm plasma, when 
the wavelength of perturbations is much larger than the Debye length. We also omit the 
negligible contribution of cosmic rays. Thus, we have 

e..= «G3 + a;,2A)^- (37) 

Substituting Equations (33)-(35) and (37) into Equation (31), we derive the following 
dispersion relation: 



^' ^ ~ (a;2 _ ^2^2) + (7^^oa;2 _ ^2^2 



^2 - _ ^2^2 ^ m pi ^ '-^^^PT^, (38) 



where we also have neglected the contribution of the displacement current. Equation (38) 
describes the streaming and thermal instability. Below, we consider some particular cases. 



10.1. Streaming Instability without Cosmic Ray Back-Reaction 

In this case, we neglect the contribution of the cosmic ray term in the dispersion 
relation defined by Equation (38). We also set all the frequencies Q to zero. Then Equation 
(38) can be written in the form 

where Cg = [7 (Teo + Tjq) /mj] ' is the plasma sound velocity, cai = oodc/oopi, pio = miUiQ, 
and jret is defined by Equation (16). Equation (39) coincides with Equation (9) in the 
paper by Riquelme & Spitkovsky (2010) (see also Bell (2005)) . 
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10.2. Streaming Instability with Cosmic Ray Back-Reaction 



Taking into account the back-reaction of cosmic rays, Equation (38) becomes 



OJ 



-k'c 



2„2 



^T^Jretk 



+ 



A'KJl.Qk 



(40) 



PiO (^2 - k'^cl) PcrO (7crOW^ " k'^C^^r) ' 

where pcro = ''^cro^^cr- Since jcro = jret, the second term on the right-hand side of Equation 
(40) is considerably larger than the first one roughly by a factor of nio/ricro ^ 1- Thus, 
the back-reaction of the streaming cosmic rays results in much more powerful instability 
than the induced background plasma streaming. It should be noted that this conclusion is 
satisfied for conditions under consideration. From Equation (40), omitting the first term on 
the right-hand side, we can find the wave number km and growth rate 6m {S = —iu) of the 
fastest growing mode 



kl 



^'^JcrO 






P^rOC' {^-Id^ - C\y 



1 



(7, 



crO scr 



>)= 



'^IcrO^scr^/ 



1/2 



1 



(41) 



and 



27^JcrO CA , -1/2 47rj2^o {icrlclcr + 4) 47crO 



PcrO^ ^scr 



7crO 



PcrOC 



\fcrO 



C^^^ ^ AJ 



.2 
'scr 



X 



'^IcrO^scr^A 



1/2 



(42) 



Let us find asymptotical expressions for km and 5m- In the case 'Jcro'^icr ^ c^, we have 



.2 _ "^^JcrO 7c/o CA .2 _ ^^JcrO 4 



2 3 ' m 2 2 



(43) 



Thus, 6m = km [icro (^scrCA) ■ In the opposite case, 7croCscr ^ ^^^ w^ obtain 



2 .,-1/2 ^^-2 

.2 ^ ^_^^^ ^ 
"^ -, «2 ^ ^ ' "^ r, /2 'crO' 



7,2 _ 4^JcrO 7crO ^2 _ ^^.^crO,,-l 



PcrO^' Cscr^'A 



PctQC 



(44) 
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The relation between 5m and km is the same as for solutions (43). From Equations (43) and 
(44), we can write expressions for k"^ and 5^, which unite both limiting cases 

7,2 _ ^^JcrO 7crO '^A ,-2 _ '^'^JcrO IcrO'^A ( AtT\ 

rn 2 / -1 9 I 2 \ ' rn o -1 9 i 2 ' V^'^ ) 

PcrOC' Cscr (TcroCLr + Ca) PcrOC' 7croCLr + ^A 

In the resonance case, 'Jcro'^lcr ~ c^j ^^ find from Equations (41) and (42) 

,2 ^ ^-^crO 1 .2 ^ ^-^crO , -1 /.^^x 

^m 2 2 ' ™- 2 'crO' V ^/ 

As we see, magnitudes given by Equation (45) in the resonance case are only twice as large 
as those in Equation (46). Thus, Equation (45) can be applied to a good accuracy for any 
relation between 7~oC^cr ^^^ '^a- 



10.3. Thermal Instability with Cosmic Ray Back-reaction 

From Equation (38), it is clear that the thermal instability is described by equation 

where C^ is given by Equation (36). If we set Ucr = 0, then Equation (47) is satisfied. In 
the case Ucr 7^ 0, the value (w^ — k^Cf) is multiplied by a small coefficient in Equation 
(38) in comparison with the second term on the right-hand side of this equation. Therefore, 
Equation (47) is kept. Thus, cosmic rays do not influence on the thermal instability under 
conditions considered in this paper. 

We set in Equation (47) Tjo = Teo = Tq. Then, this equation coincides with Equation 
(47) in the paper by Nekrasov (2011) without the inertia term. We note that in the last 
paper the perturbation of the thermal energy exchange frequency has been taken into 
account. We further take Qie = ^d = ^e- Then, we have 

S' {S' + M + /34) + Y^'^cl (27(5=^ + /3i(5 + h) = 0, (48) 
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where, as above, 5 = —ito and Cg = {2'jTo/mi) ' . The following notations are introduced: 

/3l = (7 + 1) {n^ + ^Te + ^Ti) - fine - ^ni + ^l^e, (49) 

/32 = (fix + fi^e ~ fine) fir* + {^Ti — fini) (fix + fiTe) 

+ 2 (iZ^ + ilxe ~ ^^ne ~^ fir* ~ fini) ^hy 
f33 = n^ + VtTe + fiTi + 2fi„ 

/34 = (fi^ + fire) firi + (fix + fire + fir*) fie- 

The frequencies fire,*, fine.i, fiej.je, and fi^^ are given by Equation (A7). In the general form, 
Equation (48) can be solved numerically. We note that fi^ = (7 — 1) Xeo^^/'^eo- 

We now treat Equations (48) and (49) in the limit fi^ ^ fix, fire,*, fine,*- In the short 
wavelength limit, /c^c^ ^ 5^, the dispersion relation has the form 

52 + 2fi,(5 + -fir nfie = 0, (50) 

7 ' 

where fiT,n = (fix + fire — fine + firi — fim)- Solution of Equation (50) is the following: 

5 = -^nT,n. (51) 

27 

This solution corresponds to Field's isobaric solution (Field 1965). In the long wavelength 
limit, k'^cj. <C S"^, we have equation 

S^ + 2Vt,5 + fi^fi, = 0, (52) 

where fi^ = (fix + fire + fir*)- Solution of Equation (52) is 

5 = -^fir, (53) 

which corresponds to Parker's isochoric solution (Parker 1953). 
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10.4. Thermal Instability without Cosmic Ray Back-reaction 

For ultrarelativistic cosmic rays, 7cro -^ oo, their back-reaction is absent (see Appendix 
B). In this case, Equation (38) takes the form 






6 +kC^= ^r>^-J^7J^-l7^^T^ (54) 

where 



In Equation (54), we have assumed that q^r = Qi and Ucr ~ c. We also consider that 
^jo = "^eo- In the low- frequency regime, 5"^ <C k'^c\^, Equation (54) together with Equation 
(55) is given by 



, 27(5^ + /3i(5 + /32 _ 2!^ 



^1^^ +k^^ r: :,:::, ' -^1^^. m 



If we assume that 5 ^ Wpjn^^o/''^fo) then we return to the case considered in the previous 
section. In this case, the plasma return current plays no the role. When the opposite 
condition, 6^ <C i^pj^cro/'^fo' i^ satisfied then Equation (56) takes the form 

'f;ff;f^ ^2,^'^^a. (57) 

We see from Equation (57) that in the limiting case a <C 1 (a ^ 1) the nominator 
(denominator) tends to zero. In the case a ~ 1, the dispersion relation is modified, however 
the qualitative character of the thermal instability does not change. Thus, the plasma 
return current does not influence on the thermal instability in the low-frequency regime. In 
the high-frequency regime, 5^ ^ ^^c^jj Equation (54) is the following: 

where we have assumed that 5"^ <C UpikcAincro/'i^io- Again, if 6 ^ 1 {b ^ 1), then the 
nominator (denominator) on the left-hand side of this equation tends to zero. Thus, the 



-21 - 



plasma return current has no effect on the thermal instability in these limiting cases. When 
6 ~ 1, we have solution 

d r^ Upi . (59) 

The left-hand side of Equation (58) is of the order of the unity and does not describe the 
thermal instability for the solution given by Equation (59). 



11. DISCUSSION 

In this paper, no conditions have been used for the background plasma except for 
cu^- ^ d'^/dt'^, which is usually satisfied in astrophysical settings. For cosmic rays, we 
have assumed that 7f^g d/uccrdt <^ 1 (see Section 9). This condition can be satisfied for 
moderately relativistic cosmic rays. However, it can be violated for ultrarelativistic cosmic 
rays. Using the growth rate (45) in the case for example 'Jcro'^'icr ^ '^A; ^^e last condition 
can be written in the form 



o I I fjff \ LXif 



where we have assumed that ca ~ CAi and Uccr ~ ^d- In the limit Ucr -^ c, Equation (60), 
taking into account that Tcr ^ iricrC^, can be violated for sufficiently dense cosmic rays. In 
the opposite case, 'Jcro'^'scr ^ <^i) the corresponding condition is given by 

\nioJ CAi 

We see from Equation (38) that at conditions under consideration and for one- 
dimensional perturbations along the background magnetic field cosmic rays do not infiuence 
on the thermal instability (see Sections 10.3 and 10.4). At the same time, the back-reaction 
of cosmic rays results in the aperiodic streaming instability much more powerful than that 
due to the return current of the background plasma. The maximal growth rate is achieved 
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for sufficiently cold cosmic rays and large magnetic fields such that c\ > Tcro^Lr- ^^ ^^is 
case, the growth rate is equal to 

r ^ Ucr^ -1/2 

"max ^^ ^pcr IcrO 

and the wave number is 

-1/4 
"^cr IcrO 



yCscrCA) 

Thus in particular cases such as considered here as an example, the cosmic ray back-reaction 
must be certainly invoked to study cosmic ray streaming instabilities. 

We see that the thermal instability (Equation (47)) is not influenced by the action 
of cosmic rays in the model under consideration, if we take into account the cosmic ray 
back-reaction. The multi-fluid dispersion relation includes the ion inertia and has a general 
form except for Tjo = Tgo. In the limit of fast thermal energy exchange in which Q^ is much 
larger than all other frequencies, the isobaric and isochoric growth rates have been obtained 
(Equations (51) and (53)). Ultrarelativistic cosmic rays do not experience back-reaction. In 
this case, the plasma term with the return current is contained in the dispersion relation 
(Equation (54)). However, as we have shown, this term also does not influence on the 
thermal instability. 

We note that all instabilities considered in this paper are connected with the particle 
dynamics along the background magnetic field. 

We have explored the situation in which cosmic rays drift across the background 
magnetic field. This model has been considered by Riquelme & Spitkovsky (2010) for 
the problem of the magnetic field amplification in the upstream region of the supernova 
remnant shocks. However, such a model can also be applied to the ICM where cosmic rays 
are the important ingredient (Loewenstein et al. 1991; Guo & Oh 2008; Sharma et al. 
2009, 2010). In another model, cosmic rays drift along the magnetic field. This case has 
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been investigated by Bell (2004) (see also Riquelme & Spitkovsky 2009). In both cases the 
growth rates are the same (Bell 2004, 2005; Riquelme & Spitkovsky 2010). Such a situation 
can also be encountered in the ICM. In the papers by Bell (2004, 2005) and Riquelme & 
Spitkovsky (2010), the return current of the background plasma has only been involved in 
the analytical treatments. The cosmic ray back-reaction has been included in the numerical 
analysis and found to determine the saturation of the instability (Riquelme & Spitkovsky 
2009, 2010). However, the influence of the cosmic ray back-reaction on the growth rate 
remained unknown. 

As we have obtained in this paper, the cosmic ray back-reaction drives the instability 
whose growth rate is proportional to n~Q , but not ti^q as the one due to the return 
plasma current. Therefore, this instability can produce much larger magnetic fleld 
ampliflcation in both the upstream medium of shocks and ICM. In unstable regions, an 
enhanced X-ray radiation must be observed. 



12. ASTROPHYSICAL IMPLICATIONS 

Our linear analysis of the instabilities related to the current- driven instability by 
cosmic rays is applicable to a variety of environments. Although such a type of instability 
was suggested originally for the magnetic fleld ampliflcation in the shocks by the supernova, 
we think, wherever there is a strong cosmic ray streaming, this instability may play a 
signiflcant role. For example, if the supernova driven shock propagates through a hot and 
low density medium (i.e., superb ubbles) , then the current-driven instability may exist. 
Even at larger scales, such as shocks in the ICM, we may expect this instability under 
some conditions. Most of the previous analytical studies are restricted to cosmic rays 
drifting along the ambient magnetic fleld without their back-reaction and possible thermal 
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effects. Interestingly, our analysis shows that inclusion of back-reaction will lead to a much 
stronger instability in comparison to the previous studies where this effect is neglected. 
So, we expect that the magnetic field is amplified to a larger value in the presence of the 
back-reaction of cosmic rays. This implies more confinement of cosmic rays with excited 
turbulent motions in the non-linear regime and accordingly the acceleration of cosmic rays 
to higher energies. 

In some of supernova remnants such as IC 443, SN 1006, Kepler, Tycho and etc., 
the driven shocks are propagating in their partially ionized ambient medium. This was a 
good motivation to extend CRCD instability from MHD approach to a two-fluid case, by 
considering ions and neutrals as two separate fluids where they can exchange momentum 
via collisions (e.g., Reville et al. 2007; see also Bykov & Toptygin 2005). It was shown the 
instability is getting slower rate because of collisions of ions with neutrals, in particular 
when the cosmic ray flux is not very strong. However, the back-reaction of cosmic rays has 
not been considered by Reville et al. (2007). Considering the finding that the growth rate 
is significantly enhanced in the presence of cosmic ray back-reaction in a two-fluid system 
consisting of the ions and electrons, one may naturally expect such an effect in a three-fluid 
system consisting of the ions, electrons, and neutrals. It deserves a further study, but we 
may expect that the stabilizing effect of the ion-neutral collisions will be compensated by 
the back-reaction of cosmic rays. 



13. CONCLUSION 

Using the multi-fluid approach, we have investigated streaming and thermal instabilities 
of the electron-ion plasma with homogeneous cold cosmic rays drifting across the background 
magnetic field. We have taken into account the return current of the background plasma 
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and the back-reaction of cosmic rays for one-dimensional perturbations along the magnetic 
field. It has been shown that the cosmic ray back-reaction results in a streaming instability 
having considerably larger growth rate than that due to usually treated return current of 
the background plasma. The maximal growth rates and corresponding wave numbers have 
been found. 

The thermal instability has been shown to be not subjected to the action of cosmic 
rays in the model under consideration. The dispersion relation for the thermal instability 
in the multi-fluid approach has been derived which includes the ion inertia. In the limit of 
fast thermal energy exchange between electrons and ions the isobaric and isochoric growth 
rates have been obtained. 

The results of this paper can be useful for the investigation of the electron-ion 
astrophysical objects such as galaxy clusters including the dynamics of streaming cosmic 
rays. 
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A. APPENDIX 

A.l. Perturbed Velocities of Ions and Electrons 

We put in Equation (1) v^ = Vjo + Vji, pj = pjo +Pji, E = Eq + Ei, B = Bq + Bi. For 
perturbations depending only on the 2;-coordinate, we have Vjq • V =0. Then the hnearized 
Equation (1) takes the form 

— — l-*jiH Vji X Bo, (Alj 

ot rrij rrij Ujq rrijC 

where we have used that pji = tijqTji + tijiTjq {rij = rijo + riji, Tj = Tjq + Tji) and 



introduced the notation 



F,i = ^Ei+^^v,o X Bi. (A2) 



nij rrijC 



We find from Equation (Al) the following equations for Vji^^y'- 

/^ "9^ , ,2 ^ „, r , (^Fjix 



y-Q^ + ^eij Vjlx = (^cjFjiy + -^, (A3) 



[q^ + ^c,j ^ny = -^cjF,,, + -^. 
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Applying d/dt to the z-component of Equation (Al) and using the hnearized continuity 
equation (2), we obtain 



dt"^ rrij dz'^ 



Vjl 



jlz 



rrij dzdt 



3^ I '^^ JJ-2 



dt 



(A4) 



A. 2. Perturbed Temperatures of Ions and Electrons 



Let us find now equations for the temperature perturbations. Linearized versions of 
Equations (3) and (4) for one-dimensional perturbations are given by 



Tl■^ 

DuTii = —Cii h iZjgTei, 

Tl 1 
DleTel = — Cie h fiejTji, 



(A5) 



where notations are introduced 

d 



T, 



JO 



Dli — 777 + ^Ti + ^iei Cu 

at 
d 

Die = 777 + ^Y + ^Te + ^ei, C*le = T^ 



(7-l)| + l^. 



(A6) 



dt 



eO 



-(7-i)| + a 



The frequencies Vt in Equations (A5) and (A6) are the following: 

dCi{nio,Tio) dCi{nio,Tio) 

^Ti = (7 - 1) — —^^ ' ^^ni = (7 - 1) 



UiodTi 



iO 



dCe{neo,Teo) 

"Te = (7 - 1) ;^^= ' "'^ 



Tioduio 
(7-I) — ^-^ M 



(A7) 



T^n<9n 



eO(-"ieO 



, ^u-^ 



-(7-1) 






neodTeo 

When deriving Equation (A5), we have used Equation (17) and Equations (2) and (5) 
in their linearized form. From Equation (A5), we can express temperature perturbations 
through the number density perturbations 

n ■-, n -, 

(A8) 



DTii — —DieCii VlieCie , 



DTf-l — —DiiCie VteiCii — . 
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where 

D = DuD^, - Q,,Q,i. (A9) 

To proceed further, we apply operator d/dt to Equation (A8) and use the continuity 
equation. As a result, we obtain 

D^ = D,.Cu^ + n^^cj-^^ (AlO) 



Ot OZ OZ 

These equations, we have to introduce into Equation (A4). 



A. 3. Equations for longitudinal Velocities Vn^ and v, 



elz 

Let us rewrite Equation (A4) for each component of species and use Equation (AlO). 
Then we obtain 

LuViu + L2iVelz = D^^, (All) 

ot 

T I r j^dFelz 



dt 
Here, the following notations are introduced: 

^1^ = ^7^ + — (^I'^C'i, - T,oD) --, L2^ = —n^^CieTT^, (Al2) 

ot^ rrii oz'^ rrii oz^ 

^2 1 ^2 1 a^ 

From Equation (All), we find equations for Vnz and Veiz 

LViiz = Hii, LVelz = Hei, (A13) 
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where 



L = {LuLie - L2iL2e) , (Al4) 

d 

Hil = D— [LieFiiz — L2iFelz) , 

d 

Hel = -^^ [LliFeiz — L2eFiiz) . 



A. 4. Simplification of Operators defining Vi^eiz 



Let us introduce notations 



W, = -f— + QT^ - ^ni, Vi = T^ + ^Ti, (Al5) 

ot at 



(9 n T/ — "^ 

— + iZ^ + iZj-e — iZ„e, Ve — TT: 

ot ot 



^e = 7t;I + ^X + ^re - ^ne, K = ^ + ^x + ^^e- 



Then the following operators take the form 

D = V,V, + n,iVi + (],eK, (A16) 

-l-eQ 



L = D'^ + L,D-^ + '-^^^L2 



where 



dt"^ dz^dt'^ minip dz^'' 



L, = -^ [We (Vi + n,e) + Viilei] - — [W^ (Ve + ^ei) + ^e^e] , (Al7) 

L2 = W,WeV,Ve + W,V, {We + K) ^e^ + W^eK (W^. + V,) fi,e 
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We further have 



-1 



ot J rrij \ c / c rrii 



\ ot J rUe \ c 



where notations are introduced 



df^ ITLp 



ly I tj4-Dlz, 

' c rrii 



Qe 



G2 = D^ - ^ {WeV, + W^e^e + V^ne^) ^, 
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Gg = D^^ - — 
G4 



qt 



W,Ve + W^^fie* + Ve^ie ' - {W, - Vi) il^ 

{w.-y:)^,,—-. 

nip oz^ 



dz 



dz'^' 



2' 



In Equation (A18), we have used expressions (see Equation (A2)) 



J^ilz — ^lz-\ -D 

TTli \ C 

nip \ c 



ly -Dlx I , 



l?// ' 



where Vqx = cEoy/Bo. 



(A18) 



(A19) 



(A20) 



B. APPENDIX 

B.l. Perturbed Velocity of Cosmic Rays 

For the cold, nonrelativistic, T^r <C nicrC^, cosmic rays, the hnearized Equation (6) 
takes the form 



'IcrO^E^T- + IcrO—^^cr ' ^E^]— = h F^rlH Vcrl X Bq, (B1) 



dt ^^-0^2^-- Q^ 
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where 

Feri = — ( Ei+iu,, X Bi ) (B2) 

and Ucr is directed along the y-axis. We have used that jcri = 7cro^cr ■ ^cri/c^, where 
7cro = (1 ~ '"cr/'^^) • Equations (Bl) and (B2) do not include Veo in Equation (13). From 
Equation (Bl), we find the following equations for Vcrix,y'- 

IcrOl^+^ccr ^crlx = UJccrFcrly + IcrO^K—, (B3) 



crO pi,2 ccr j ^crlx ^ccr-^ crly ~r IcrO o, ) 

4 0^ . 2 \ _ IP I U-tcriy 
'ycrO'^'^'^ccr I ^crly — —^^ccr^crlx + IcrO — ^T ■ 

The z-component of Equation (Bl) is given by 



7crO a^ = ^ h Fcrlz- (B4) 

From Equation (7) in the linear approximation, we find 

/ Tlcrl 2 '^crVcrly \ f-nr\ 

Perl = PcrOi^cr 7crO o • (B5) 

\ncrO c J 

Applying to Equation (B4) operator d/dt, substituting Equation (B5), and using the 
continuity equation for cosmic rays, we obtain 

Lcr^crlz = Hcrl- (B6) 

Here, 

Lev = 7crO^ - '^scrg^^ (B7) 

2 2 ^cr (y ^crly , ^^ crlz 



-fJcrl — f^scr IcrO 9 ci cu. ~r 



where Cscr = {Pcro^ cr / iT^'cr'ncro) is the sound speed of cosmic rays. We note that the first 
term on the right-hand side in the definition of i^cri in Equation (B7) is connected with the 
perturbation of the Lorentz factor. 



